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A dielectric function including collisional correlations is derived by linearizing the self consistent 
Vlasov equation with a Fokker-Planck collision integral. The calculation yields the same type of 
dielectric function as in the standard-theory of Schottky-noise in storage rings. This dielectric 
function is compared with the Mermin-dielectric function derived from a kinetic equation with 
relaxation time approximation. We observe that these functions are identical, however the Mermin- 
DF is computationally advantageous. The limits of both dielectric functions are given and the sum 
rules are proven. We apply these dielectric functions for typical storage ring plasmas and calculate 
the stopping power and the plasmon excitation spectrum. 



I. INTRODUCTION 

During the last ten years experiments with ions ( p + 
up to U 92+ ) in storage rings gained importance in the 
field of spectroscopy and plasma physics. The storaged 
and cooled ion beams have a high luminosity for recom- 
bination experiments and inertial confined fusion inves- 
tigations. In particular it is of basic interest to study the 
transition between the weak and strong coupled plasma 
or even the transition to a crystalline state of a cooled 
ion beam Q. The most important prerequisite for ob- 
taining dense states is strong electron and laser cooling. 
The electron cooling force can be described as stopping 
power acting on an ion beam in an electron plasma 
Another plasma phenomena in dense beams are collective 
excitations (plasmons, shear modes) which are detectable 
by the Schottky noise ||. All items - the pair distribu- 
tion function of a state, the stopping power and the shape 
of the collective excitations are related to the dielectric 
function e(q, m). 

Within the linear response theory the polarizability 
H(q, lo) (and altogether the dielectric function e(q,Lu)) is 
defined by the variation of particle density Sn(q,Lu) in 
reaction to an external field 5U ex t(q,uj) via 



U(q,u) 



su ext 



(1) 



The connection to the dielectric function (DF) is given 
by 



e(g» = l + V c U(q,w) 



(2) 



The captured ions in storage rings are moving in front 
of a background of the confining fields ensuring approxi- 
mately the charge neutrality in the system. This nearly 
neutral system of ions interacting via the Coulomb po- 
tential Vc immersed in a homogeneous background of op- 
posite charge is usually called an one-component-plasma 
(OCP). 



An unpleasant problem is the temperature anisotropy. 
The longitudinal temperature (Tj|) differs from the 
transversal (Xj_) (referring to the beam-axis) because 
only the longitudinal direction is cooled directly. Even 
taking into account relaxation processes between the de- 
grees of freedom and possible transversal cooling the tem- 
perature difference maintains. 

In this paper we focus on experiments done by the 
Heidelberg TSR group with 9 Bc + and 12 C 6+ ions cooled 
by an electron beam jjjy]. The Be-ions can be cooled 
further down to few mK by applying laser cooling. 

For estimations about the ideality and collision- 
numbers we employ the longitudinal temperature Xji . 
The density of the ion beam can be calculated using the 
current j of the beam profile (the diameter) measured 
by the beam profile monitors x bpm and the value of the 
betatron function on this position (3b pm 



jQPbpm 
2ttZx 2 bpm Rvq 



(3) 



Here denotes vq the ion beam velocity, R is the ring radius 
and Z is the charge state of the ions. The transversal 
tune Q amounts to 2.8. 



Parameter 


Be+ 


C b+ 


e 


n 


2.3 • 10 la m-' 5 


1.6 • 10 i3 m -3 


2.9 ■ lO^m™ 3 


T|| 


6K 


8000K 


3.5K 


r 


0.13 


0.0031 


0.23 


A 


1.7 • 10's -1 


4.0 • 10 s s _i 


3.4 • lO^s" 1 




2.1 ■ W°s' L 


9.2 ■ lO^s -1 


3.0 • lO'V 1 



TABLE I. Parameters for several experiments in the TSR. 
The explanations are given in the text. Data are from [111 
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The essential parameter for characterizing the coupling 
is the nonideality or plasma parameter 



r = 



e 2 f Ann Y 
Ane Q k B T \~) 



(4) 



which is the ratio of potential and thermal energy. Fur- 
ther essential quantities are the plasma-frequency 




J P i 



and the inverse screening length 



(5) 



The validity of these sum rules is an essential statement 
about the quality and reliability of the dielectric funct ion. 
The outline of the paper is as follows: In Chapter [I A 



we give a short rederivation of Mermin DF and th e DF 
from a Fok kcr-P lanck equation is given in Chapter 



II B 



In Chapter [II C| we co mpare both DF's a nd th e sum rules 
are proven in Chapter [I E . In Chapter III A we present 
the application to the stopping power and in Chapter 
IIIB we calculate the Schottky spectra. 



II. DIELECTRIC FUNCTIONS WITH 
CORRELATIONS 



(6) 



Another important parameter characterizing the plasma 
is the collision frequency, or friction coefficient A which 
is the inverse relaxation time 



A = - = na(v th )vth 

T 



(7) 



cr{y) = An 



(- 



2e 2 



\12e k B T 



A(v). 



Here is <r(v t h) the cross-section at thermal velocity v 2 h — 
2k B T '/to for ion-ion collisions, A(u) denotes the Coulomb 
logarithm e.g. in Brooks-Herring approximation. 

The collisions between the ions play an essential role 
in this storage-ring-plasmas being responsible for effects 
like intra-beam scattering (IBS) la]. That means an ex- 
pansion of the (ion) beam due to ion-ion collisions has 
to be suppressed by electron cooling. Consequently for a 
sufficient description a dielectric function e(<f, oj) includ- 
ing these collisions should be considered || . In this paper 
two practical ways will be shown in chapter 

An essential property of every dielectric function is the 
fulfillment of the sum rules. The most strongest are the 
longitudinal frequency sum rule 



— ^-Ime 1 (q, uj)duj = 1 



and the conductivity sum rule 
2uj 



Ime(<f, uj)Auj = — 1 



o % 



moreover the compressibility sum rule 

'°° 2 1. 



lim 

9^0 



o 



7TUO U> 



— Ime (q,u))d(J = 1 



and the perfect screening sum rule 



lim 



2 1 

— 2 — Ime(q, uj)duj 



= -1. 



(8) 



(9) 



(10) 



(11) 



A. Mermin-dielectric function 

For calculating a dielectric function including collisions 
between the particles with mass to, N. D. Mermin 
suggested a particle number conserving dielectric func- 
tion. We shortly sketch his derivation for the classical 
case starting with a kinetic equation in relaxation time 
approximation 

ot or or to ov 



f(f,v,t)-f (v) 



(12) 



This kinetic equation describes the development of a par- 
ticle distribution function f(r, v, t) consisting of an equi- 
librium part fo(v) and a non-equilibrium part Sf(f,v,t) 



f(f,v,t) = f (v) + Sf(r,v,t) 



(13) 



The mean field U(f, t) is composed of an external part 
and a part arising from the induced particle density dn 



U(q~iv) = 5U ext + VSn(q, uj) 



(14) 



One gets the induced particle density Sn by linearization 
of ( p^ ) and integrating the solution of 5f over the velocity 
v. After Fourier transformation t — > u> and r — > q the 
following polarization function Tl(q, to) is obtained 



5n(q,cj) = 



$f(q,v,uj)d 3 v 

no(£w + i) 



-5U ext {q,Lu) (15) 



l-V(0IIo(&w + i) 
with the RPA or Lindhard polarization function 

1 #/(«) 

771 OV "V/ 



n (<7» = / d a v 



vq 



(16) 



The RPA dielectric function 

.2 



UJ + i/r 
i in classical limit reads 



e(<7,w) 



1 



hi' 



1 — 2x r e 



2k B T q 



(17) 
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and fulfills all sum rules (p|)-(|ll|). Shifting the frequency 
into the complex plane according to [lfj] one gets the relax- 
ation dielectric function. This expression does not fulfill 
the limit of static screening and has a non Drude-like 
high frequency behavior which leads to a violation of the 
sum rules (f)-©. 

In jjj was suggested a more sophisticated dielectric 
function by considering the relaxation ansatz 

|/(f, v, t) + v§-J(f, v, t) + JLfft ^ t) = 

Ot or or m Ov 

f(r,v,t) - f (r,v,t) 

T 

with respect to a local equilibrium distribution function 

-~ fj, + 6/j,(r,t)~ 



fo(r,v,t) = cxp 



mv 
2k B T 



k B T 



(19) 



instead of the global distribution fo(v) in (jl2|). In the 
simplest case one can specify the local distribution by a 
small fluctuation in the chemical potential <5/i related to 
the density-fluctuation Sn. 

The Mermin-dielectric function is derived by solving 
( |l8| ) using an expansion of the local equilibrium distribu- 
tion function in powers of Sfi 

f (q, v, lo) = f (v) - qdvf °S V) 6n(q, w) . (20) 
mqv 

Sfi is determined by the particle number conservation 
ojSn(q,ui) = J vqSf(q,v,uj)d 3 v. leading to @ 



5n{q,uj) 

n(<f,o) 



(21) 



Finally one obtains from ( |18[ ) and (|20|) for the polariza- 
tion function 



n M (g,w) 



n (g,w 



1 - 



1 



1 - 



n (<f, u) 



(22) 



l-iwr \ n o (q,0) J 

Instead of (|l^) we arrive at a density variation 
n M (£w+i) 



l-V(q)Hu(q,uJ + ^) 



5U ext {q,uj) 



(23) 



and with (^|) the Mermin-dielectric function finally has 
the shape 

- | = 1 + V utI V W t> . (24) 

lut e(q, 0) - 1 



e M [q,v + 



Here denotes e(g, u> + ^) the dielectric function ( |l7| ) in re- 
laxation time approximation. It is easy to see, that in the 
limit r — > oo the Mermin-dielectric function reproduces 
the RPA-dielectric function (|l7|). 



B. The Vlasov-Fokker-Planck Equation 

Now we examine another kinetic equation - the Vlasov 
equation with the Fokker-Planck collision integral which 
has been used to predict the Schottky noise of an ion 
beam [|| 

d d d U(r,t) d 

Tr:f(r,v,t) + v — f(f,v,t) + — — — f(r,v,t) 

at or or m ov 



=4 (xlr + ^ 



(25) 



The application of the Fokker-Planck collision term is 
valid for weak collisions (it means low g-values) because 
it represents an expansion of the collision integral in mo- 
mentum space. With the collision integral of the Fokker- 
Planck equation one includes the fluctuations of the dis- 
tribution function due to collisions. It describes the bal- 
ance between dynamical friction A (w/(r, v, t)) holding 
the velocity-distribution sharply near zero-velocity and 
the diffusion D-^f(f, v, t) flattening the velocity distri- 
bution. The coefficients A and D in the Fokker-Planck 
equation are related by the Einstein relation 



D 

7 



k B T 



(26) 



As already mentioned above the friction coefficient A is 
equal to the inverse relaxation time. Obviously the drift 
coefficient Xv is linear in the velocity as long as the dif- 
fusion coefficient D is a constant. The Fokker-Planck 
collision term ensures the particle conservation. Due to 
the Einstein relation we have a proper balance between 
friction and diffusion. So we expect that similar physics 
is included as the Mermin extension of the simple relax- 
ation time approximation in the last paragraph. 

We solve this Fokker-Planck equation again within the 
linear response. A sketch of the derivation can be found 
in appendix A with the result for the dielectric function 



e VFP (q,Lo) = 1 + — 



k B T q 2 
m X 



ito 



k B T 



1 i , k BT 2 .W 

1,1 + - l y m \^ 



(27) 



and iFi denotes the confluent hypergeometric function. 
This dielectric function has been given in [|| and is valid 
for an isotropic plasma in three dimensions. 
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co[s 1 ] 

FIG. 1. Comparison of the real parts of the RPA- (circles), Mermin- (squares) and VFP-DF (diamonds). We have chosen 
temperatures of 0.6 ■ • ■ 60K available in the longitudinal direction of an ion beam and wave numbers below and above the inverse 
Debye length k (^). The particle density is n = 2.3 ■ 10 13 m -3 of single charged Beryllium ions. The real parts of Mermin- and 
the VFP-DF are identical. 
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FIG. 2. Comparison of the imaginary parts of the RPA- (circles), Mermin- (squares) and VFP-DF (diamonds). The same 
density and temperatures are chosen as in the Fig. pi The imaginary parts of Mermin- and the VFP-DF are identical. 
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C. Comparison of both dielectric function 

Up to now we have used different kinetic equations 
leading to two different dielectric functions. It is of great 
interest how these functions are related to each other and 
whether these dielectric functions are valid in the storage 
ring plasma's realm of temperature, density and friction 
coefficient. In Fig.|j] and Fig. || are plotted both dielectric 
functions in dependence on the frequency for several wave 
numbers and temperatures. 

We see that within the numerical accuracy of the pic- 
ture no difference is visible between the dielectric func- 
tion of a Fokker Planck collision integral and the Mermin 
dielectric function resulting from the conserving relax- 
ation time approximation. 

Let us inspect now some special limits. Both dielectric 
function fulfills the static limit (ui — > 0) 



e(g,0) = 1 + — 
q 2 



(28) 



for all A in accordance with the classical Debye-Hiickel 
result for static screening. In the long-wavelength limit 
q — ► one gets the Drude formula for both dielectric 
functions 



lim e(q,u) = l pl . . 

q^O u(u + iX) 

co — ► oo 



(29) 



For A — ► this formula reproduces the RPA behavior. 
In the limit of strong friction A — * oo we get in concur 
with H and || also (|2^). The long- wavelength and the 
strong friction limits are identical. 

For low temperatures there are differences between the 
RPA-dielectric function and the other correlated dielec- 
tric functions. The real parts start in the static limit at 
the same value as the RPA-dielectric function but drops 
down much more earlier (in Fig.^ one sees 4 magnitudes 
at one tenth of the inverse Debye length and 2 magni- 
tudes at ten times the inverse Debye length). There are 
no zeros in the real part. Accordingly, the imaginary 
part is shifted in the same fashion. It is one magnitude 
broader than the RPA imaginary part and has only two 
thirds of its height. For temperatures higher than 50K 
the RPA-dielectric function and the Mermin- and VFP 
dielectric functions become identical. 



D. Sum rules 

The most interesting question is whether the dielectric 
function fulfills the sum rules (|)-(0). Due to (|§) all 
presented dielectric functions lead to Hee(q, u>) oc lu^ 2 for 
large u>. Since poles due to the relaxation time occur only 
in the lower half plane we have 



dui , K _ ' . = 

(tJ + IT]) 



(30) 



from which we see that the dielectric functions fulfill the 
Kramers-Kronig relations 



Ree -1 ^) -1=V 



2 Ime 1 (q, to) duj 



2tt 



(31) 



where the V denotes here the Cauchy principle value. 

From (|3l| ) we get with in the static limit just the 
compressibility sum rule (lid). The longitudinal f-sum 
rule (||) follows as well from (pl]) . To see this we observe 
that due to time reversibility e(q,—u>) = e*(q,u>) holds 
and we can write 

OG 

lim Ree _1 (g, w) = I + lim [ dwlme _1 (<7, u>). (32) 

ui — >oo uj — >oo UJ^TT J 



Using ( p9| ) we obtain just the f-sum rule ( ph. 

Since the same Kramers-Kronig relation (|3l|) holds also 
for e instead of e _1 we see that the corresponding free sum 
rules (||) and ([ll]) are also fulfilled. 

This completes the proof that both correlated dielec- 
tric functions fulfill the sum rules. We can state therefore 
that the dielectric function is properly valid in the inter- 
esting scope and can be used to describe the phenom- 
ena in cold and dilute storage ring plasmas. Since the 
Mermin dielectric function is computational much easier 
to handle, in fact no more circumstances than the Lind- 
hard dielectric function, we will use the Mermin dielectric 
function further on. 



III. APPLICATION TO STORAGE RING 
PLASMAS 

We continue now to apply the correlated dielectric 
function derived in the last paragraph to typical storage 
ring plasmas. Two important quantities we like to dis- 
cuss here. First the stopping power of ions in an electron 
plasma and second the occurring plasmon excitations. 



A. Stopping power 

The stopping power, i.e. the energy transfer of a parti- 
cle to a plasma is given in terms of the dielectric function 

by [| 



dE a 



d*q 



at h j (2^)3 huj Vaa{q)2 Ime ~%M. 

(33) 

Here rib denotes the Bose function and V aa the Coulomb 
potential of the particle a. We observe that the sum 
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about different plasma species is condensed in the dielec- 
tric function. It is noteworthy to remark that this result 
is valid for any arbitrary degeneracy. The derivation pre- 
sented in M shows that the result ( |33| ) is more general 
valid than derived earlier |Io|-|l2|] . Higher order correla- 
tions like vertex corrections can be incorporated in the 
dielectric function, such that ( |33| ) remains valid Q . This 
fact is important for dense solid state plasmas which are 
used recently for stopping experiments, where the result 
( |33| ) is applicable as well. A more explicit form can be 
given by carrying out the angular integration [q = Tik] 



v(t)k+4 



dE a 2e 2 a 1 fdk 



8t 7T£ 



1 f dk f 
—-^ — / duj uj iib{u) Ime _1 (?ifc, oj). 
v{t) J k J 



(34) 



Neglecting the quantum effects in ( j34| ) which represent 
an internal ionic cut-off due to the thermal De Broglie 
wavelength we get the standard result of dielectric theory 



dE a _ 2e 2 a 1 fdk 
dt ireo v{t) 



v(t)k 



JjT J duuj line- 1 {Tik, u) (35) 
o o 

from which all known special cases can be derived Jl , 
among them the well known Bethe formula. We use (B4) 
where no artificial cut-off is needed further on. 

In figure || we have plotted the stopping power of 9 Be + 
calculated with the Mermin and the Lindhard dielectric 
function. We observe that for a weakly coupled stor- 
age ring plasma with temperature of 6 K and a den- 
sity of 2.3 10 13 m -3 which corresponds to a nonideality 
of r = 0.13 almost no differences are observed between 
the Mermin and Lindhard result. For higher coupling by 
lower temperature of 1 K corresponding to L = 0.77 we 
see that the Mermin stopping power become smaller than 
the Lindhard result. Since the friction is dependent on 
the squared of density but only on temperature via the 
Coulomb logarithm we find a stronger dependence on the 
density. This is illustrated in the figures [3|||. We see that 
with increasing density the deviations between Mermin 
and Lindhard results become appreciable. 
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FIG. 3. The stopping power of 9 Be + ions in an elec- 
tronic plasma versus ion energy. The classical Lindhard 
result (thick lines) are compared with the Mermin result 
(thin lines) for two different temperatures. The plasma 
parameters are V — 0.77 (solid lines) and F = 0.13 
(dashed lines) respectively. 
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FIG. 4. The stopping power of 9 Be + ions in an elec- 
tronic plasma versus ion energy. The classical Lindhard 
result (thick line) is compared with the Mermin result 
(thin line) and the T-matrix result (dashed line) of bi- 
nary collisions (p6|). The plasma parameter is F = 1.65. 
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9 Be* energy [eV] 

FIG. 5. Same situation as in Fig. ^, but lower density 
T = 0.77. 



ical results for the transport cross section for a plasma 
with charge Z = 1. In figures [| - ^| we compare the re- 
sults for the dielectric theory of stopping power with and 
without collisional contributions with the pure two- par- 
ticle collision result of ([56|). We see that the two-particle 
collision expression is significantly smaller than the di- 
electric theory. For very strong coupling in figure ^ we 
see even a vanishing contribution of the latter one indi- 
cating that the two-particle collisions do not contribute 
any more but the energy transfer is mainly caused by 
collective spreading. In Fig.R we represent the reduced 
energy loss [Ai = e 2 /(127re T)] 

Xivth 1 dE 
ksT v dx 

versus the coupling parameter ZT 3 / 2 . 



(38) 
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FIG. 6. Same as in Fig. but density is now 4 
magnitudes higher (r = 16.5) 

So far we have generalized the dielectric theory of stop- 
ping power by the inclusion of collisions. It is instructive 
to compare now the results directly with the stopping 
power in binary collision approximation. In the follow- 
ing expression for the stopping power was derived from 
the Boltzmann equation within T-matrix approximation 



dE 
dt 



(«)=E 



n b v t 



a cosh a — (1 



p 2 (l 



m a kBT 



-)sinh a 



'U+T^) 2 



(36) 



with the thermal velocity v 2 — 2kBT/rrib, the abbrevi- 
ation a = jr^p{l + and the quantum mechanical 
transport cross section 



n-' {l> ) /dfi(l-COBfl)^. 



(37) 
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In fl3fl a fit formula is given which subsumed the numer- 



FIG. 7. The normalized friction coefficient (energy loss) 
for RPA(solid line), T-matrix result (dashed line) and Mermin 
DF result (dot-dashed line). The filled circles are simulation 
results which reproduces experimental data |lE| , [L6| ] . 

The dependence of the normalized energy loss from the 
coupling parameter is weaker in the Mermin case than in 
the RPA case but distinct from the numerical simula- 
tions. Nevertheless the involving of collisions modifies 
the stopping power in the right direction. The best de- 
scription is still given by the T-matrix result (pq). 



B. Plasmons 

The ion beam current j(t) is a fluctuating quantity 
due to its granular (ionic) structure. Detecting the 
mirror charge on the vacuum chamber of the ring and 
Fourier transforming (frequency analyzing) one obtains 
the Schottky-signal. It is primarily used for analyzing 
the beam's velocity distribution and hence the longitu- 
dinal temperature, but also important for measuring the 
particle number or the revolution frequency of the beam. 
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It is related to the dynamical structure factor S(q, uS) by 
the equation |L7[ 



\J(q,oj)\ 2 ) ~S(q,u>) 



(39) 



where the brackets indicate the thermal averaging. The 
well-known fluctuation-dissipation theorem |l8f| connects 
the imaginary part of the response function lme(q, u>) and 
the dynamical structure factor 



S{q,u) = 



k B T 
LuV C {q) 



Ime 1 (q,uj) 



with the Coulomb potential Vc- In dense ion beams (e.g. 
the C 6+ LIR experiment) one observes a double peaked 
Schottky spectrum. The two peaks commonly identified 
as plasma waves propagating in two directions around 
the storage ring. A frequency analysis of this beam cur- 
rent shows the propagating waves clearly as peaks in the 
spectra which are theoretically well described in H|). 
The identity of Mermin and VFP-dielectric function we 
use now to compute the the Schottky noise much easier 
within the Mermin DF. 

For numerical calculation one has to modify the plasma 
frequency which differs from the one of an isotropic 
plasma. For a plasma in a conducting tube we have 



Cu 2 



pi 



27ri?Qmeo7 



In^ + i 

r B 2 



(40) 



Here denotes N the particles number, tq the radius of 
the beam chamber, rg the radius of the beam propa- 
gating with the velocity vo = jc (speed of light), 2nRo 
is the circumference of the ring. M is the number of 
plasma waves fitting in the ring, the so-called harmonic 
wavenumber. 

The equation for the plasma frequency follows straight- 
forward from inserting the beam pipe's impedance Z via 
Ohm's law U = e 2 voZ/Ro in the external field of the 
Vlasov equation. Usually Z is given by 



Z = 



M 
vol 2 



In 



(41) 



We assume for simplicity the non-relativistic case which 
is valid in the TSR experiments ( 9 Be+: 0.04c, 12 C 6+ : 
0.15c) with 7 « 1. 

The wavelength q is not a continuous variable but as- 
sumes only discrete values q = M /2ttRq . The fraction 
K 2 /q 2 can now expressed by 



q 2 5u 2 

with the thermal frequency 

2k B T 
R 2 m 



(42) 



Scj z = 



Inserting this parameters into the VFP dielectric function 
one obtains 



e V FP {M,lo) 



6u> 



1F1 



1,1 



M 2 Slo 2 
~llX 2 ~ 



iuj M 2 5lo 2 
T' 2A 2 



(44) 



This is the well known standard permittivity in the 
Schottky noise theory. Since the Mermin DF and the 
VFP-DF are identical we use practically now the easier 
Mermin DF. We modify equation (24) according to the 
parameters (|40|)-(p|) 



e M (M,Lu + i\) = 1 + 



(1+|) (e(M,a; + a)-l) 
| iX e(M, w + iX) - 1 



u e{M, 0) - 1 
with @ for e(M, ui + iX) and 



e M (M,0) = 1 



,3l 

'Slu 2 



) (45) 



(46) 



(47) 



In the next step we insert a relaxation time considering 
the anisotropy in the thermal velocities Vth 



A = 4ir(ZeY 



e 2 m 2 v th ^v 2 hll 



A. 



(48) 



Here A again denotes the Coulomb logarithm. We have 
used the modified Mermin-dielectric function for calcu- 
lating the plasmonic excitation for a 12 C 6+ beam at an 
energy of 73.3MeV and a revolution frequency of vq/Rq — 
617kHz (M — 5). In || the expression Ime(M, ui)/uj is 
compared with the Schottky measurement. We see that 
the modified Mermin-dielectric function fits the Schottky 
spectra satisfactorily. 
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FIG. 8. The Schottky spectra of a dense carbo- 
nium beam and the corresponding theoretical prediction. 
( 43 ) (n = 8.3 ■ 10 13 m" 3 , Tu = 11000K). Data are taken from 111. 
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IV. SUMMARY 

In this paper we have described two dielectric func- 
tions including collisions. After numerical inspection we 
have shown that the DF obtained from VFP equation 
is identical with the Mermin DF. Modifying the dielec- 
tric function for storage ring purposes (modified plasma 
frequency, discrete wavenumbers/harmonics) we have de- 
rived the standard dielectric function for Schottky noise 
prediction for a three dimensional plasma beam. Because 
of the identity of both dielectric functions one can use 
a Mermin dielectric function for Schottky noise descrip- 
tion, too. The second goal was a better description of 
the stopping power acting on an ion beam in the cooler's 
electron gas. Here we like to state, that including the 
collisions leads to a lower friction force than in the RPA- 
predictions but obvious overestimates the friction force 
compared with simulated and experimental results. 

Further improvement has to be done to consider strong 
coupling effects in the relaxation time by using a match- 
ing Coulomb logarithm. Efforts have also to be done 
to include the magnetic field in the cooler and the 
anisotropic Maxwell distributions in the dielectric func- 
tion. 



One can identifying the functions H„ as Hermite- 

Polynoms. The lowest eigenfunction for 

is the Maxwell-distribution. The eigenvalues are 



A 



-\{n x + n y + n z + 2k 2 ) 



(A5) 



In the next step we insert for the right side the eigen- 
values and eigenfunctions of the inhomogeneous Fokker- 
Planck operator 

»q$At)-™m,t)yp- (A6) 

at m ou 

= A £ fW-^w. 

With fin(K) = — - — 
A 

and expand the distribution function Jg/o and the distor- 
tion 6f(u, k,t) in a sum of eigenfunctions of the Fokker- 
Planck operator 

- d - °° 

k-7^fo(u) = k ^ a ni n v ntll>n inv n, («, k) 



Sf(u,k,t)= ^2 CnxnynA 



(A7) 



APPENDIX A: SOLUTION OF THE 
VLASOV-FOKKER-PLANCK EQUATION 

We following the the main ideas of pfj|| . At first we 
introduce reduced variables 



r <f k BT 



(Al) 



2k B T 

After linearization the VFP-equation ( p5|) we arrive at 

^-6f(k, u, t) + 2i\uk 5f(k, u, t) - —kU(k, 

at m ou 



X— (- — 

du V 2 du 



(A2) 



For the Fokker-Planck term including the second term of 
the left side (the so-called inhomogeneous Fokker-Planck 
operator) one can consider the eigenvalue and eigenfunc- 
tions in three dimensions 



-^Xuk + X-^--^ 



(k, u) (A3) 



The eigenfunctions take the form after a coordinate 
transformation z = u + 2ik 



( m 



V njn y \n z \2 n * +n « +n * \2irk B T / 
e k2 e~^-^ 2 R nx [z x }U ny [z y ]U nz [z z ] . (A4) 



Performing a Fourier transformation t — > to we arrive at 
a solution for the coefficients 



^kanUik,^) 



y —iuj — \{n x + n y + n z + 2k 2 ) 
Remembering (|15|) we obtain 



(A8) 



8ll(k,u)= ^ c n x n y n, J VVi^n^U (A9) 



ii'xii-yii'z Vol 

= U(k,ui)HyFp(k,ui). 
Hence we get for the polarizability 

nvFP(fc,cH^-A^exp[fc 2 ] J2 i 2 ^ +n " +n ^ 

r|| + n x ] r_[| + n y ] n z ] T [| + n x + n y ] 

l(n x + n y )\(n x + n y + n z )\ 
(n x + n y + n z + 2k 2 )(V2k) 2 ^+ n y +n ^ 

x " w JTT^ ■ ( A1 °) 

— iu — X(n x + n y + n z + lk z ) 

Using the relations 

r [i + n x ] r [i + n y ] r [\ + n z ] r [± + n x + n y ] 



E 



n x ,n v ,n* 



+ fly + 71*)! 



E- 



(All) 



and 







exp[— x] 



n\(n + n) 
we arrive finally at 

n V Fp(fc,^) = 

1 



iFi[1,1 + k,:b] (A12) 



(A13) 



2k B T 



2Xk 2 



— 1F1 



l,l + 2fc 2 -i-;2k 
A 



The dielectric function related to this polarizability 



was discussed in chapter [IB 
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